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Orbiting detectors will be able to observe showers initiated by neutrinos 
penetrating the Earth and interacting close to their exit point. There is 
a correlation between the impact parameter of the incident neutrino and 
OO ! its energy. We study the development of upward going, neutrino induced 

showers in the atmosphere. 
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: 1 Introduction 

> ■ 

^ . The study of extraterrestrial, ultrahigh energy neutrinos is important from at least two 

points of view. First, they probably carry astrophysical information about point sources 
which are optically thick at (almost) all wavelengths, for instance, an active galactic 
nucleus (AGN). Second, incoming neutrinos of laboratory energies of the order of 10 18 eV 
(or even higher) are expected to be emitted by the highest energy sources. In a collision 
with a nucleon, the CMS energy is, therefore, of the order of 40 TeV, far exceeding the 
CMS energies available in accelerator generated neutrino beams. Although such events 
are expected to be rare ones, they provide a unique window on physics beyond the 
Standard Model. For both reasons it is important to refine the observation of ultrahigh 
energy (UHE) neutrinos. The main purpose of this paper is to sketch a novel way of 
observing such interactions. 

The basic principle is a very simple one. Neutrino-nucleon interactions have a cross 
section which is a well known, monotonically increasing function of the incident neutrino 
energy; for a recent calculation, up to neutrino energies of 10 12 GeV, cf. ref. [[L|. Conse- 
quently, on the average, a neutrino penetrating the Earth will interact after a well defined 
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distance beneath the surface. If the interaction occurs reasonably close to the exit point 
of the neutrino, the ensuing shower develops largely in the atmosphere and it can be 
detected by means of an orbiting detector, such as OWL (Orbiting Wide angle Light 
collector) or AIR WATCH. Such detectors are suitable for a study of the longitudinal 
development of showers. Thus, a shower developing upwards gives a unique signature of 
a neutrino interaction. Further, one obtains some information about the primary energy 
of the incident neutrino: the Earth acts as an energy filter. Indeed, given the impact 
parameter of the incident neutrino, its (energy dependent) interaction mean free path 
has to be approximately equal to the length of the chord of the trajectory inside the 
Earth. If it is considerably shorter (high energies), the interaction and the development 
of the shower takes place largely inside the Earth: no atmospheric shower is observed. 
Conversely, at low energies, the neutrino penetrates the Earth without interacting. Once 
it leaves the Earth, its chance of interacting within the atmosphere is rather slim. 

This paper is organized as follows. In Section ^| we summarize the calculation used 
in order to obtain the shower development. The results are given in Section |3|. Finally, 
in Section f| the results are discussedQ. A preliminary account of these results was given 



We use one dimensional cascade theory in Approximation A; due to the fact that we are 
mainly interested in ultrahigh energies, this should be adequate in order to obtain the 
first estimate of the behavior of the showers induced by neutrino interactions. 

The incoming neutrino interacts with a quark in the target nucleus and generates 
a hadronic jet. In the present work, we assume that the hadrons generated share the 
primary energy eqally; hence we start with showers induced by a single hadron. 

In a hadronic interaction, most of the secondaries produced are pions. Fast nucleons 
are relatively rare, with the possible exception of the leading nucleon in a nucleon-nucleus 
interaction. At this point, we neglect the production of particles other than pions. 

We assume that pions of either charge are produced in equal numbers. Since the 
average multiplicity in an interaction is large, this is a permissible simplification. We 
neglect photoproduction of pions: hence, neutral pions act as the source of the electro- 
magnetic component without the latter affecting the hadronic one. As a consequence, 
the equations governing the hadronic component are decoupled from the ones governing 
the electromagnetic one. 

We assume that Feynman scaling is valid and parametrize the single particle inclusive 
distributions by means of a simple expression, 



Here, z = E/E' is the Feynman parameter, E' being the primary energy. The quantities 
2 Throughout this paper we use natural units, i.e. ti = c = 1 
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K,a,f3 are to be determined from the normalization conditions, 



dE 



da 



(2) 




and 



ldEE% = E 
J dE 



(3) 



and from fits to the experimental data. An infrared (IR) cutoff, z Q has been intro- 
duced. Equation (0) follows from the definition of a single particle inclusive cross section, 



Due to the fact that we treat all light hadrons together, some compromises are neces- 
sary: in fact, the parameter {3 governs the leading particle behavior which is different for 
the various secondaries. Nevertheless, the available data || suggest that taking a = 
and P = 5 gives a reasonable fit to all data. (By choosing a = one gets an inclusive 
distribution oc \jz for small values of the Feynman parameter. This is consistent with 
the fact that at low momenta the inclusive distribution is dominated by soft gluon emis- 
sion.) The infrared cutoff has been chosen at a fixed Feynman parameter rather than at 
a fixed energy. This is consistent with the fact that in the following, we use an energy 
independent average multiplicity, < n >= 50 — again a compromise. By doing so, one 
underestimates the number of soft secondaries. However, we found that these approx- 
imations have virtually no effect on the integral distibutions with a threshold energy 
E^^lOOGeV. Likewise, the shower development is rather insensitive to the value of /3: 
even a distribution F{z) oc 0(1 — z) gives reasonable showers. Likewise, the integral 
distributions at the threshold energies indicated above are quite insensitive to the choice 
of the IR cutoff. For this reason, we put z = whenever this choice does not give rise 
to IR divergences. 

We are now ready to write down the transport equations for the hadronic component 
of the shower. There are two distinct equations. The first one is for the stable hadrons: 
nucleons and charged pions; at the relevant energies the charged pions may be considered 
stable. The second equation characterizes the behavior of neutral pions: in that case, 
decay is essential^. We denote the differential distribution of charged hadrons and neutral 
pions by h and p, respectively. The depth of the medium is measured in units of the 
hadronic interaction mean free path, r = t/X, where t is the target thickness measured 
in g/cm 2 . The approximate transport equations read as follows. 



In writing down equation (f|) we neglected the interaction of neutral pions. The loss 
term contains merely their decay. The factor - serves to convert the real space decay 

3 One recalls that in air at ground level, the decay and interaction mean free paths of a neutral pion 
become equal at an energy of about 10 18 eV. 
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h + 2/3 f^F(§)h(E') 
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length, D, to the units used in writing down the transport equation. The decay length 
is, of course, given by the expression: 



TE 

where T is the lifetime in the rest frame and mo is the mass of the neutral pion. The 
density of the medium, p, is expressed as a function of r. The mass of the neutral pion 
is denoted by m . The factors 2/3 and 1/3 reflect the fact that approximately two out 
of three hadrons produced are charged pions and one out of three is a neutral pion. 

In a similar fashion, we can write down the equation governing the electromagnetic 
component. To the accuracy desired for the purposes of the present calculation, the 
(common) inclusive distribution of Bremsstrahlung and pair production can be approx- 
imated by a step function, F(E/E') = 0(1 — E/E'). (This gives, of course, a pure 1/E 
inclusive spectrum for both the charged leptons and the Bremsstrahlung photons, a rea- 
sonably accurate approximation.) In the same approximation, we put the cross sections 
of Bremsstrahlung and pair production equal to each other (i.e. 7/9 ~ 1). We denote 
the number of particles interacting only electromagnetically by I, i.e. 

n e + + n e - + n 7 ~ I. 

The latter approximate equality is valid after a few radiation lengths; initially, the num- 
ber of charged leptons is slightly overestimated. 

We write down the the evolution equation for I by measuring distances in units of 
the radiation length, X . Assuming again that all neutral pions decay instantaneously 
and that the energy of the 7T° is shared equally between the decay photons, we have with 
o = t/X : 

^ = §p (2B, + / §W - E)l (£', „) - I (E, a) (5) 

The evolution equations given by and (|5]) form a hierarchical structure. Due to the 
(very reasonable) approximation made in writing down these equations, namely that 
neutral pions decay instantly and do not interact, one has to solve only the first one of 
equations @; the evolution equation for neutral pions can be solved by quadrature and 
equation (^) is solved in terms of a Green function. 

In fact, the evolution equation for neutral pions is solved by the expression: 

p = J dr'e- {T - T ' )x/D Q(r - r')g(r'), (6) 

where g(r) is the gain term as one can read it off from the second of equations (|J). Due 
to the fact that X/D ^> 1, one can pull out g(r') from under the integral in eq. @ 
at t' = t. Upon substituting this into the evolution equation of leptons, eq. (|5J), one 
realizes that in the rapid decay approximation, that equation also scales as the hadronic 
evolution equations are supposed to do. Hence we conclude that Feynman scaling holds 
for a mixed hadron - lepton cascade if: 



4 



• The energies are sufficiently high so that effects of the rest masses of the particles 
are negligibly small, 

• the energies are sufficiently low so that the reinteraction of decaying neutral par- 
ticles is negligible together with the QCD loop corrections containing logarithms 
of the energy. (If the decay length does not cancel, there is a scaling violation due 
to the Lorentz factor present in the expression of the decay length.) 

In order to find the retarded Green function of eq. (|5|), one uses Mellin transform 
techniques. The resultis: 

G(q,r-r')=exp(r-r')Q(r-r')Q(l-q)(^^j ' 

x/ 1 (2 3 / 2 (r-r') 1/2 (ln|g|) 1 /2). (7) 

In eq. (|7|), q stands for the scaling variable, q = E/E' and J is a modified Bessel function 
of order one, cf. 

To summarize: assuming the validity of Feynman scaling and the rapid decay of 
neutral pions, one has to solve one evolution equation only, viz. the one governing the 
evolution of charged hadrons. Other components of the mixed hadronic and leptonic 
cascade are obtained by quadratures, which is rapidly accomplished by using standard 
numerical integration techniques. 

The solution of the first of the equations @ is best accomplished by iteration. In 
fact, one notices that by converting that equation into a Volterra integral equation, 
an iterative procedure rapidly converges. Roughly speaking, the number of iterations 
needed is about the number of collision mean free paths at which the distribution is to 
be computed. (This has been known for a very long time: Bhaba's method of "succesive 
collisions" for obtaining the evolution of an electromagnetic cascade is, in essence, based 
on this observation, see, e.g. 0.) Equations for the integral distributions can be derived 
in a straight forward manner; the calculation is particularly easy in the Feynman scaling 
limit. 

3 Results 

We calculated the integral distributions of electrons and positrons for a number of hadron 
energies and nadir angles. In each case, we assumed that the incident neutrino interacts 
in the Earth, just below its exit point. (In practical terms, this means that the interaction 
has to take place within a few hadronic interaction mfp below the surface.) For each nadir 
angle, the target thickness penetrated by a neutrino before interaction was determined on 
the basis of the results given in ref. 0. The profile of the electromagnetic component of 
the shower is shown in the following Figures. In each calculation, the integral spectrum of 
the hadrons was cut off at E = lOOGeV, corresponding to a CMS energy of y^s « 14GeV. 
Below such energies, particle production becomes insignificant; moreover, the particles 
produced tend to decay into soft muons and gammas instead of contributing to the 
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cascade. Likewise, the integral distribution of the electromagnetic component was cut 
off at E c = lOOMeV, roughly corresponding to the critical energy in air. In estimating 
the primary neutrino energies, we assumed an average inelasticity < y >= 0.5 and an 
average hadronic multiplicity, < n >= 50 in the first interaction. Thus, on the average, 
Ehadron = E u /100. The variation of the nadir angle around the horizontal direction 
corresponds to a slight variation of the target thickness around one neutrino mfp. 

The following figures display the electron-positron component generated by a single 
hadron. Due to the linearity of the evolution equations, the curves displayed here can be 
easily scaled to accommodate other inelasticities and/or neutrino primary energies. The 
shower profile is presented in "real space". For purposes of converting target depth to 
length an exponential atmosphere was used, with scale height, h = 16.8km and ground 
level density, p — 1-3 x 10~ 3 g/cm 3 . The notation used in the Figures is summarized in 
Table 1. 



E h (GeV) 


Nadir angles 




diamonds 


triangles 


asterisks 


10 7 


80° 


85° 


88° 


10 s 


80° 


84° 


87° 


10 9 


80° 


84° 


88° 



Table 1. Summary of hadron energies and symbols used for showers developing at 
different nadir angles. The average single hadron energy, Eh is assumed to be 1% of the 

incident neutrino energy, see text. 
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Fig.l: Integral distribution of electrons and positrons, E > lOOMeV. E h = 10 7 GeV 
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Fig.2: Same, E h = 10 8 GeV 
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Fig.3: Same, E h = 10 9 GeV 



4 Discussion 

The figures displayed in the preceding section reveal a characteristic feature of the show- 
ers induced by a netrino interaction near the exit point from the Earth. Those showers 
are "upside-down" - they start in the dense portion of the atmosphere and evolve up- 
wards. As a consequence, such showers reach their maximum rather soon, typically at 
distances of the order of 5km. Thereafter, they have an unusually long tail, i.e. a very 
slowly decreasing particle number as the shower evolves. This is an effect which is easily 
understood: the "upside-down" shower starts near the surface of the Earth and proceeds 
upwards. As far as the hadronic (charged pions and nucleons) and the purely leptonic 
(7, e^) components are concerned, the patterns of a "normal" or of an "upside-down" 
shower are the same if plotted against the target depth. However, the distribution of neu- 
tral pions (the main source of photons responsible for initiating the leptonic component) 
is asymmetric: upward and downward going beams of 7r°-s can be easily distinguished 
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from each other. As a consequence, any orbiting detector which is able to follow the 
longitudinal development of a shower is, a fortiori able to pick up upward going showers. 

It remains to be seen whether sufficiently powerful extraterrestrial neutrino sources 
exist so that one can collect a significant number of neutrino events of the type described 
here. Besides the obvious astrophysical interest of such events, one would be able to look 
for signatures of "new physics" , assuming that it affects the longitudinal development of 
the neutrino induced showers. For instance, it was proposed that at some higher energy 
scale, neutrinos begin to interact more strongly than it would be expected on the basis of 
the standard electroweak theory see refs. || |7|(]. Should this be the case, the correlation 
between the impact parameter of the incident neutrino and the development of the 
shower would change. Assuming a sufficiently large statistical sample of "anomalous" 
events, one could deduce the onset of a new physical phenomenon^. 
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4 It was pointed out that, taken literally, such schemes violate unitarity, see This is due to the 
fact that in refs. (|, Q the onset of "new physics" was modelled by adding a piece to the cross section 
which is proportional to a step function. As a consequence, the real part of the corresponding amplitude 
develops a logarithmic singularity. Almost any smooth function approximating a step function is free 
of such a singularity and can be reconciled with unitarity. 

5 Of course, fluctuations can affect the results; this question is being studied. However, at sufficiently 
high energies, the development of a shower gives information on the energy of the shower on an event 
by event basis as well, see, e.g. 
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